Applicability of Fermi golden rule and possibility of low-field runaway transport in 

nitrides 
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In order to justify applicability of the standard approach of perturbation theory for the description of 
transport phenomena in wide-band polar semiconductors with strong electron-phonon interactions, 
we have compared dependences of energy losses to the lattice on the electron drift velocity obtained 
for different materials in the frameworks of (a) a perturbative approach based on calculation of 
the scattering rates from Fermi's golden rule and (b) a non-perturbative approach based on the 
path-integral formalism of Thornber and Feynman. Our results reveal that despite strong electron- 
phonon coupling in GaN and AIN such that intercoUision times become of the order of the period 
of phonon oscillation, standard perturbative treatment can still be applied successfully for this type 
of material. Our findings also indicate possibility for unique long-distance runaway transport in 
nitrides which may occur at the pre-threshold electric fields. Polaron ground state energy and 
effective masses are calculated for GaN and AIN as well as for GaAs and AI2O3. 

PACS Numbers: 72.80.Ey, 72.10.-d, 72.10.Di, 63.20.Kr 
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I. INTRODUCTION 



One of the main factors determining electron trans- 
port characteristics in polar semiconductors is scatter- 
ing of the electrons by polar optical phonons. For rel- 
atively weak electron-phonon interactions, when scat- 
tering events can be considered as independent, use 
of Fermi's golden rule for the calculation of energy- 
dependent frequencies of electron transitions provides an 
adequate description of experimentally obtained velocity- 
field curves. Upon increasing of the interaction, however, 
the polaronic effects induced by autolocalization of an 
electron by the inertial part of the crystal polarization be- 
come more prominent and they determine the character 
of the scattering. Intensification of the electron-phonon 
interaction leads eventually to a situation where the av- 
erage intercoUision time becomes less than the duration 
of a collision. Such a strong coupling, therefore, requires 
proper account for the quantum interference effects and 
makes the problem of electron drift essentially nonlinear. 
This complicates dramatically the theoretical treatment 
of carrier scattering and field-dependent transport since, 
for the given case, the standard perturbation technique 
is not applicable. 

To ensure energy conservation for the short-time per- 
turbations, the inverse scattering rate, t, must be large 
enough to satisfy the inequality r ^ 7i/Af , where A£ is 
the electron transition energy. This criterion, however, 
does not allow one to conclude whether or not a wide 
class of materials in which the transition energy may 
become of the order of the linewidth can be described 
successfully in the framework of standard perturbation 
approaches. To such a class of materials belong, in par- 
ticular, nitrides of Al and Ga. These semiconductors 



have been investigated recently quite intensively due to 
a number of their unique properties that can be utilized 
in the current state-of-the-art semiconductor technology. 
However, most attempts to describe scattering processes 
in nitrides have been undertaken assuming the validity 
of Fermi golden rule. Moreover, despite the possibility 
of growing the nitrides of group III in zincblende-like 
structures, their crystal structure at ambient conditions 
is wurtzite-like. For wurtzites, it is generally necessary 
to account for optical anisotropy when considering the 
carrier-optical-phonon interactions. Since optical phonon 
spectra in wurtzites are far more complicated than those 
of cubic crystals, the majority of theoretical results have 
been obtained by ignoring the features of the phonons in 
optically anisotropic media. Recently, a formalism has 
been developed^ for evaluation of the scattering rates in 
bulk wurtzite-like semiconductors and heterostructures 
by taking into account peculiarities of the phonon spectra 
obtained in the, framework of macroscopic dielectric con- 
tinuum modelfl Due to the anisotropy-induced complex- 
ity of the problem, Ref. |l| took advantage of the pertur- 
bation theory. However, as indicated by the previously- 
discussed considerations, the validity of such an approach 
requires independent confirmation. 

The present paper demonstrates applicability of 
Fermi's golden rule for describing adequately the 
electron-longitudinal-phonon interaction in polar materi- 
als by comparing the field- velocity dependences obtained 
in the frameworks of (a) the perturbation theory and (b) 
the non-perturbative path-integral approach of Thorn- 
ber and FeynmanEI (TF). A supplemental, but very im- 
portant, result obtained from the present investigation is 
the discovery of the possibility for unique low-field long- 
distance runaway transport in materials characterized by 
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strong electron-phonon interactions. Examples of such 
materials are the nitrides of Al and Ga. 



II. MODEL 

The most systematic and self-consistent approach for 
evaluating the long-range polaron ground-state energy 
G, effective mass toq, and carrier energy dissipation for 
both strong- and weak-coupling limits of electron- ph^ 
interaction has been developed by Feynman et al 
In Ref. ^, the problem of electron drift in a parabolic 
band under steady state conditions is considered quan- 
tum mechanicallyQ assuming that all the energy losses are 
due to interaction of electrons with polar optical modes. 
Taking advantage of the Frohlich's polaron model, the 
authors used the path-integral method to eliminate the 
lattice coordinates from the momentum balance equa- 
tion and obtained an expression for the magnitude of the 
electric field E that is required to maintain a particular 
magnitude of electron velocity V at arbitrary tempera- 
ture and interaction strength characterized by the cou- 
pling constant 
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Herein, e is the elementary charge. Evaluation of a re- 
quires four parameters: electron effective mass m*; fre- 
quency of the longitudinal phonon Q; static dielectric 
constant eo; and high frequency dielectric constant eoc- 
All these parameters can be measured experimentally 
and they are the only external parameters required for 
calculation of energy loss ner unit distance eE versus V. 

Thornber and FeynmannI have calculated the depen- 
dences of eE{V) for three coupling constants (a = 
3, 5, 7) over a wide range of reciprocal temperatures 
f3 — h^l/ {ksT), where fcs is the Boltzmann constant 
and T is the temperature of the lattice. The general 
result of these calculations can be summarized briefly 
as follows. For each particular a, eEiV) has a maxi- 
mum at some threshold value Vth- For /3 > 1, location 
of this maximum becomes independent of temperature. 
For V < Vth, sE is an increasing function of velocity. 
This interval of velocities corresponds to a stable situa- 
tion when energy loss to the lattice due to the absorption 
and emission of optical phonons can be compensated by 
the energy gained by the electron from the applied field 
in such a way that at the given E, a small deviation AV 
of the velocity from its steady state value Vg creates a 
force, e[E{Vs) — E{Vs ± AF)], which stabilizes the veloc- 
ity at Vg. When the external field approaches the value 
Eth = E{Vth), the dependence tends to saturate since 
the magnitude of the energy loss due to interaction with 
optical phonons is finite. The case E > Eth was ex- 
cluded from consideration because no steady state con- 
ditions can be reached for such fields and electron would 
accelerate infinitely. The theory, however, predicts the 



existence of solutions for V > Vth- In this region, eE 
is a decreasing function of the velocity which leads to 
an unstable steady state situation. For this case, any 
deviation of the velocity from Vg would lead to either 
deceleration of the electron to velocity V < Vth which 
is stable at the given field, or a gradually increasing ac- 
celeration if AV leads to an increase in velocity. It is 
essential, that for V > Vth, the dependence eE{V) can 
be interpreted as a time-dependent momentum loss in 
the absence of the external field. This loss would co- 
incide with the rate of electron momentum loss if the 
criterion dV/dt ^ V^/ (duration of the collision) is satis- 
fied. For /3 > 1, eE{V > Vth) also becomes independent 
on temperature. 

In order to simplify the comparison, and taking into 
account that the strongest electron-polar-optical-phonon 
scattering is due to emission of the longitudinal optical 
(LO) phonons, we will consider the case when /3 « 4. Due 
to high energy of LO phonons in the nitrides, such value 
of P would correspond to room temperature in these ma- 
terials. For GaAs, which we take as a reference point 
in our investigation, (3 — 4: would correspond to lattice 
temperature of order of 104 K. Since f3GaAs is slightly 
higher than 1 at the room temperature, the result ob- 
tained for maximum energy loss per unit distance can be 
compared to the experimental velocity-field dependences 
(see, for example, Ref. ^ and citations therein). Indeed, 
at some threshold field Eth, the dependence V{E) has a 
maximum caused by transitions of the carriers to an up- 
per valley with a higher effective mass. In terms of the 
TF model, these transitions would start to occur when 
the energy supply from the external field would exceed 
the maximum loss to the lattice; i.e., at E > Eth- 

Thus, if the average kinetic energy of electrons ob- 
tained in the framework of TF model at Eth does not 
exceed the energy of bands separation and the effects of 
the increased effective mass due to non-parabolicity of 
r band can be neglected, the value of the argument at 
maximum of V{E) dependence has to correlate with the 
extremum of the eE{V). 

Under the assumptions made above, we use the sim- 
plest model for estimation of energy loss to the lattice in 
the framework of the perturbation theory. In this model, 
we calculate the dependence of scattering rate 1 /r due to 
the emission of LO phonons on a single electron kinetic 
energy £ using the Fermi's golden rule. For an optically 
isotropic material, the scattering rate is 
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where Nq is the phonon occupation number. Then, we 
assume that a carrier with velocity V = ^2£jm* loses 
the energy ?iO to the lattice in a distance tV . 
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III. RESULTS AND DISCUSSION 

A. Limiting cases 

The energy loss per unit distance vs. the electron ve- 
locity calculated for GaAs at room temperature in the 
framework of the model which uses Fermi's golden rule 
is shown on Fig. |^ by the thin solid line. It is important, 
that as anticipated, the maximum of the dependence is 
in a good agreement with the experimental data which 
give the maximum of V{E) at E fa 3.4 — 3.9 kV/cm. One 
should note, that since calculations are made in a one- 
electron approximation and for single parabolic band, 
they overestimate the velocity obtained at the maxi- 
mum. The experimentally measured values at maximum 
of V{Ep reflect the averaging of the velocities over bands 
with different effective masses as well as effects of non- 
parabolicity. 
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FIG. 1. The energy loss per unit distance vs. electron ve- 
locity in GaAs. Perturbative model; room temperature - thin 
solid line; T=104 K - thick solid line. TF model: T=104 K 
- upper dashed line; T « 20 K - dash-dotted line. Corrected 
TF model: T=104 K - lower dashed line. Polaron parameters 
are given in the inset. 

The curve calculated in the simplest model at 104 K 
(/3 = 4) is given by the thick solid line. Surprisingly, the 
eE(V) dependence computed for /3 = 4 in the framework 
of the TF model (upper dashed line) exhibits a maxi- 
mum located at somewhat higher fields. We assume that 
such a discrepancy occurs because in a weak-coupling 
limit the zero-order distribution for the electrons in this 
model reduces to a drifted quasi-Maxwellian. An essen- 
tial requirement for such a distribution to be valid in 



the given case is the presence of high electron concen- 
tration and strong electron-electron interactions, which 
provide randomization of the direction jsf electron mo- 
mentum between the scattering events.l3 It is unlikely, 
however, that such a randomization can be achieved for 
scattering with emission of polar optical phonons. In- 
deed, in the weak-coupling limit the motion of the car- 
rier in the near-threshold fields becomes essentially one- 
dimensional. Conservation of energy and momenta, valid 
in the weak-coupling limit, require the angle between the 
carrier and scattered phonon momenta to be no more 
than a,rccos{hil / £) . Due to this condition, emission of 
an LO phonon at threshold fields causes deviation of the 
direction of electron momentum from the direction of ap- 
plied field by no more than « 20 degrees. Additional fo- 
cusing of the electron momentum in the direction of the 
electric field comes from the inverse dependence of thja 
interaction matrix element on the phonon wavevectorO 
and leads to overestimation jof energy losses when using 
the Maxwellian distribution.t2l 

In order to resolve this discrepancy, one can suggest - 
in analogy to the classical case - that at the same electron 
temperature, reduction of the dimension would corre- 
spond to a reduction in the average carrier energy. Since 
the deviation of the direction of the electron momentum 
relative to the direction of the electric field is small but 
finite, for the weak-coupling limit we have reduced the 
energy scale in TF model by a factor of two in order to 
match the maxima. This corresponds to decreasing the 
velocity and energy losses by factors 2^/^ and 2^/^, re- 
spectively. The corrected curve eE{V) is shown by lower 
dashed line. 

The differences in the shapes of the curves obtained 
in the simplest model and the TF model occur for the 
following reasons. The former model considers only the 
emission mechanism, whereas the latter model also takes 
into account absorption. As shown in the figure by dash- 
dotted line calculated for the uncorrected case of TF 
model for /3 w 21, elimination of phonons absorption by 
reducing the lattice temperature yields the same slope 
of eE{V < Vth) as in the simple model. As expected, 
for this temperature interval, the temperature decrease 
does not affect the shape of the curve at and beyond the 
maximum. For V > Vth, the discrepancy in eE{V) be- 
tween the models appears to be due to the relatively high 
value of dV/dt in the unstable region. In the framework 
of the TF model, one can estimate this value by eE/mo, 
where mo can be obtained as uiq « vrCv^ juP", v and 
LO are the parameters of TF model. We have computed 
these parameters from minimization of the free energy at 
zero temperaturcfl Our estimations show that assuming 
the duration of the collision to be equal to r, the value 
of the derivative would be much less than V Jt only for 
V ^ 10^ cm/sec. Thus, for V > Vth the simple model 
cannot be used for GaAs under the conditions when an 
external field is applied. 

Additional confirmation of the idea that the correc- 
tions required for the application of the TF model in the 
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weak-coupling limit are induced by focusing of the car- 
rier momentum comes from the fact, that in the case of 
strong coupling - when the directions of momenta are 
randomized due to strong electron-phonon interaction 
discussed previously - the model correctly explains the 
experimentally-obtained results. In Fig. y we depict the 
energy losses calculated for AI2O3. The maximum on the 
dependence obtained in the TF model (dashed curve) is 
in the excellent agreement with the experimentallv_ipb- 
tained maximum losses in this material, 0.03 eV/AO 
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where lolzi ^z, ^Lt, and ujt are the characteristic fre- 
quencies of the Al(LO), Al(TO), El(LO), and El(TO) 
modes, respectively. The phonon frequency as a func- 
tion of 9 can be obtained from thr: dispersion relation for 
the extraordinary bulk phonons.EI The dependence a{d) 
calculated for GaN is shown on Fig. ^. 
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FIG. 2. The energy loss per unit distance vs. electron ve- 
locity in AI2O3 at room temperature. Perturbative model: 
solid line; TF model: dashed line; Polaron parameters are 
given in the inset. 

As expected, the maximum losses calculated in the per- 
turbative model (solid line) are less by an order of mag- 
nitude. Of course, since for the given material a > 1, 
the simple model is not valid and we have presented here 
both dependences simply as a means of estimating the 
possible error which can be induced by a perturbative 
treatment and to demonstrate that no energy scale re- 
duction can fit these dependences. 



B. Intermediate case: nitrides 



Due to the optical anisotropy inherent to wurtzites, 
the coupling parameter of polaron theory a becomes de- 
pendent on the angle 9 between phonon wavevector and 
the optical axis. Assuming that = ^ define this 
dependence as 
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FIG. 3. Dependence of the polaron coupling parameter on 
the angle between the phonon wavevector and the optical axis 
in GaN. 

To obtain the eE{V) dependence in the TF model we 
have used a — 0.46. This value corresponds to the enerajs 
of LO phonon calculated for GaN in the cubic phase.t2l 
The energy of the LO phonon in cubic AIN is taken to 
be 113 meV which is - as for the GaN case - between the 
energies of Al and El LO modes in wurtzite phase.H The 
scattering rates are calculated according to the formal- 
ism developed in Ref. |l|. Comparison of the dependences 
obtained for the nitrides is given in Fig. ^ Again, the 
maxima obtained in the simple model correlate very well 
with the threshold electric fields of ^(fi)-, dependences 
computed in the Monte Carlo techniquetS for a three- 
valley model for the conduction band: 140 kV/cm for 
GaN and 450 kV/cm for AIN. Note that in order to match 
the maxima one needs to use the same reduction of the 
energy scale when calculating the energy losses in the TF 
model as for the case of GaAs. Additionally, one can see 
that the shapes of the dependences for V > Vth are al- 
most the same. This agreement between the simple and 
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the TF models is due to the extremely short duration of 
the collisions which can be estimated roughly as the in- 
verted scattering rate, r ~ lO"^** sec. The increase of the 
polaron effective mass cannot compensate the increase in 
the energy loss near the threshold value and, therefore, 
cannot reduce dV/dt. Nevertheless, due to frequent col- 
lisions, the criterion dV/dt ^ V/t is satisfied for the 
nitrides even at ^ — > V,t. 



Indeed, let us consider a GaN sample in an external elec- 
tric field of 100 kV/cm. As shown in the figure, the two 
steady-state solutions for electron energy would corre- 
spond to such a situation. One of the solutions lies in 
the region V < Vth and, consequently, it reflects a stable 
solution with respect to the electron energy fluctuations. 
Another one falls into the unstable area, V > Vth- Sup- 
pose an electron is injected into the sample with an en- 
ergy somewhere in between the energy which corresponds 
to the second, unstable, solution and the threshold en- 
ergy. 




FIG. 4. The energy loss per unit distance vs. electron ve- 
locity in GaN and AIN calculated at the room temperature. 
Perturbative model: AIN - upper solid line, GaN - lower solid 
line; TF model: AIN - dash-dotted lines, GaN - dashed lines. 
Thin and thick broken curves are the uncorrected and cor- 
rected dependencies, respectively. Polaron parameters are 
given in the insets. 

The comparison made here allowed us to conclude that 
in materials for which tD, > 1, application of Fermi 
golden rule to explain transport phenomena is as good 
as in materials traditionally handled with the perturba- 
tion theory, i.e., in the materials for which much stronger 
criterion, rfi 3> 1, is satisfied. 



C. Low- field runaway transport 

It is important to emphasize another result of our in- 
vestigation. Figure || represents the dependencies of en- 
ergy losses on electron kinetic energy obtained in the 
framework of the TF model at room temperature for the 
three materials considered: GaAs, GaN and AIN. On the 
figure, vertical arrows indicate the energy of the closest 
upper valley in the corresponding conduction band. The 
picture shows clearly the possibility to achieve unique 
pre - threshold - field runaway transport in the nitrides. 
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FIG. 5. The energy loss per unit distance vs. electron ki- 
netic energy. Vertical arrows indicate the energy of the closest 
upper valley in the corresponding conduction band. T — 300 
K. 

In this case, since the energy losses to the lattice would 
exceed the energy gain from the external field, the elec- 
tron would decelerate until the stable solution at given 
field would be reached. If, however, the energy of the 
injected electron would just slightly exceed the value of 
high-energy steady-state solution, the electron would ac- 
celerate moving downwards on the unstable branch of the 
dependence until it gains enough energy, £* , to appear 
in the nearest upper valley of the conduction band. The 
value of acceleration would gradually increase due to the 
increasing difference between the force caused by the ex- 
ternal field, F = 100 keV/cm, and the force caused by 
energy losses to the lattice represented by the eE{£) de- 
pendence. In order to estimate the minimum runaway 
length Lr , we have assumed that the energy of the in- 
jected electron is 0.7 eV and that the maximum force, 
F — eE{£*) , is constantly applied to the carrier. Under 
this assumption, we obtained L^'^^ > 220 nm. 

It is easy to see that because the intervalley separa- 
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tion in AIN is smaller, the effect in this material is ex- 
pected to be not as strong as in the previous case. As- 
suming injection energy 0.31 eV and applying field 300 
kV/cm, we get L^'^ > 39 nm. Our results also show 
that the previously-discussed runaway transport cannot 
be achieved in GaAs due to the small intervalley gap and 
broad peak on the eE{£) dependence. 

The results presented in this paper were obtained ig- 
noring the non-parabolicity effects. These effects, how- 
ever, would not change our main findings qualitatively. 
In order to improve accuracy of the expected quanti- 
ties, further investigation is required. It is interesting 
to note, that the abstract possibility of the low-field run- 
away transport was mentioned initially by Thornber and 
FeynmanH Nitrides of Ga and Al promise to be materials 
where such transport could be actually realized. 



IV. SUMMARY 



In present paper, we have compared the energy losses 
to the lattice calculated in different polar semiconduc- 
tors within the frameworks of both non-perturbative and 
perturbative approaches. Our results reveal that applica- 
tion of Fermi golden rule for calculation of the scattering 
rates in nitrides, where tD, w 1.5, is as appropriate as ap- 
plication of this standard perturbative treatment in the 
materials for which the well-known criterion ril ^ 1 is 
satisfied. This finding dramatically simplifies the analysis 
of transport phenomena in the wide-band polar semicon- 
ductors with intermediate magnitude of polaron coupling 
factor, a < 1. 

Applying the non-perturbative path-integral approach 
of Thornber and Feynman to evaluation of field depen- 
dent electron energy dissipation in AIN, GaN and GaAs, 
we have found that pre-threshold low-field runaway elec- 
tron transport can be realized in the nitrides. The con- 
ditions for such a transport can be formulated as follows: 
(a) the energy of the injected carrier should exceed the 
energy which corresponds to the solution of the momen- 
tum balance equation located on the unstable branch of 
eE(V); and (b) the separation between the energy of in- 
jected carrier and the energy of the bottom of an upper 
valley must be high enough to provide a finite value of 
runaway length. It must be at least a few times higher 
than the energy of the polar optical phonon. 
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